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Abstract 

Our last three papers |S [3 El provide an occasion to make some brief notes on 
ideology of waves in plasmas and to rehabilitate Vlasov prescription to calculate rele- 
vant logarithmically divergent integrals in the principal value sense. In this approach 
asymptotical solutions of plasma oscillations obtained with the method of Laplace 
transform must be in the form of a sum of wave exponents with amplitudes being se- 
lected according to self-consistent boundary physical conditions. Landau damping is 
absent in this case by definition. Boundary electrical field together with conditions of 
absence of unphysical backward waves (for boundary problem in half-infinite plasmas) 
and of kinematical waves define single-valued dependence of boundary distribution 
function on electron velocity v in the case of transversal waves and on the surface 
break of the normal electrical field in the case of longitudinal oscillations. We have 
proposed physically more justified modified iteration procedure of collisional damp- 
ing calculation and demonstrated some results of damping decrements calculations 
in a low-collision electron-ion plasma. We have pointed out a possibility of natural 
selection of boundary conditions, which allows to eliminate simultaneously both back- 
ward and kinematical waves as for transversal as well as for longitudinal oscillations. 
Dispersion smearing of both longitudinal and transversal high-frequency waves, for 



which the smearing decrement 5x is proportional to Au>/{u>yuj'^ — ujf^) (where Au! is 
the frequency interval of exciting boundary field), and non-Maxwellian background 
distribution function might be the main causes of waves amplitude damping in colli- 
sionless plasmas imitating Landau damping. 

PACS numbers: 52.25 Dg; 52.35 Fp. 
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After the basic paper of L. Landau in 1946 which stated that even in collisionless 
Maxwellian plasmas electron waves are damping ("Landau damping"), there appeared an 
enormous amount of works, including all textbooks, that use it as a conception. In practice 
this scientific trend can be considered as being completed without any hope to carry therein 
something essentially new. 

Landau had considered the simplest case of plane half-infinite geometry with the before- 
hand given form of a solution exp{ikx) f {t) , where for determination of the function f{t) he 
had used Laplace transformation, and had postulated to calculate an appearing in Laplace 
image fp logarithmically divergent integral not in the sense of the principal value (as it 
had been proposed in 1937 by Vlasov P| for the analogous integral appearing at solving his 
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"Vlasov equations" ) but according to Landau rule of by-passing poles of (if ^-integrand in 
the complex- valued plane Vx {vx is an x-component of the electron velocity, x is a direction 
of the wave traveling). 

However it is well known (see, for instance review [3] or Eq. (1) in |4j) that asymptotical 
Landau solution in the form of a solitary plane damping wave exp{ikx—iujt) with u = ujQ—i6 
does not satisfy Vlasov equations neither at these nor at any other real (if considered as a 
limit at tending imaginary parts of uj, k to zero, see, for instance, or complex values 
a; and k. In spite of some paradoxes and hardly to be explained consequences of Landau 
theory (see, for instance, [Hj) it was supposed that Landau solution corresponds to reahty 
and appears to be true, so that Vlasov equations must be corrected by additional terms, 
which are determined by Landau rule of by-passing poles (see, for instance, |S1IE])- 

Nevertheless one can formulate a problem of finding asymptotical solutions of Vlasov 
equations and also equations of propagating transversal electromagnetic waves (with the 
same difficulty of logarithmically divergent integrals) without these additions, which appear 
to be artificial. So, the paradoxes of Landau solution are thought to be generated by 
Landau representation of logarithmically divergent integrals and, in the main, by the single- 
exponent form (suggested by Landau as well as by Vlasov, however with no paradoxes in 
the last case ) of an asymptotical solution. 

Such asymptotical solutions have been written in our last papers jH HI |H] proceeding 
from the following propositions: 

1. The logarithmically divergent integrals appearing at Laplace transformation have to 
be calculated in Vlasov prescription of the principal value sense. 

2. For evaluation of principal values of integrals we use a series of successful approxima- 
tions of the type 
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After this the obtained asymptotical solutions with the collisional damping take the 
form of a sum of exponents of the type exp(±2A;x±2C(jt + X(5) and exp(±2fcx=Fict;t — a;5), 
(5 > 0, corresponding to double poles pi, p2- 

3. The given boundary condition E{0,t) (for instance, in purely boundary problem) 
defines also, in a single-valued way, the boundary (and initial) conditions for the per- 
turbed distribution function fi{v,0,t), fi{v,x,0) according to linear integral equa- 
tions following from the condition of absence of kinematical waves 0. It guaranties 
the physically justified proportionality of the boundary and initial distribution func- 
tion to the boundary electric field. 

4. Selecting free boundary parameter Fp-^, being a Laplace transform of dE{x,t)/dx\x=o 
one can exclude (cancel) unphysical (for the case of half-infinite plasmas) backward 
waves contained in the solution in the form of residua exponential sum. These waves 
are divergent at x — > cxd in the case of a low-collision plasma. 

5. In the case of longitudinal waves at the given boundary field E{0, t) the boundary self- 
consistent function fi{v, 0, t) is determined analogously from the condition of absence 
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of kinematical waves. But for elimination of the unphysical backward waves in this 
case one must inevitably assume a break in the normal constituent of electrical field 
£^(0, t) at plasma surface (see jH] for a way to calculation its value). 

6. In general case the solution can not be represented as a single exponent even if the 
boundary condition E{0,t) has such a form. At the boundary field frequency u the 
asymptotic solution must contain at least two simultaneously existing forward wave 
exponents exp{±iujt =p ikx), corresponding to the frequencies +uj and —u. 

7. Amplitudes of all modes with different kn^u) corresponding to the given exciting field 
£^(0, t) are strongly correlated. Cancellation of kinematical and then backward waves 
must be achieved for all modes in the general asymptotical solution. 

8. Solutions of the wave equations must be real- valued whenever the boundary and initial 
conditions are such ones. 

Only after proceeding these procedures one can calculate strong relations between ampli- 
tudes of the different forward waves including, in the general case, both electron and also 
ion and hybrid ion-electron branches. 

This allows to avoid all paradoxes at calculation of E{x,t) and paradoxical tangling of 
the distribution function in with generating some strange electrical field. 

When all these conditions being fulfilled one can easily construct, using the method 
of two-dimensional Laplace transform usual solutions for the poles of image Ep^p^ , corre- 
spondingly, the long before known dispersion relations k{uj) 01IZIIH]. Thus, for the case 
of collisionless plasma mystic Landau damping and kinematical waves are really absent 
by definition. The solution appears as a sum of exponents with amplitudes selected in 
accordance with boundary and physical conditions. 

We also have developed the more physically justified modification of iteration process 
in the presence of collision terms of kinetic equation (low-collision electron- ion plasmas). 
By this method we have obtained damping decrements as for electron longitudinal waves 
as well as for two branches (low- and high-frequency ones) of transversal electromagnetic 
waves with the unusual^ decrement for the low-frequency branch. The coUisional dissipative 
absorption in the high-frequency branch aX u > is proportional to 1/ y^l — uj'i/uj'^ and 
grows to infinity sX uj ^ ujl [H]. 

In this connection we think that Van Kampen waves (see [3]) are called to compensate 
the erroneous consequences of Landau theory. 

It is highly believed that Landau damping is detected and verified experimentally. In this 
respect we should note the relatively small number of such verifications and the necessary 
extreme delicateness of similar experiments. The theoretically absent Landau damping in 
fact might be imitated by a series of secondary effects. These ones can be: 

(i) the difference of the background distribution from Maxwellian one due to electron 
collisions with and recombinations on the walls of discharge tube; 

(ii) the presence of longitudinal magnetic field and cyclotron motion and transversal dif- 
fusion to walls; 

(iii) effects of the method of plasma oscillations exciting; 

(iv) effects of reflecting from walls; 
^nonlinear in concentrations of charged particles. 
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(v) effects of the base electric field supporting discharge; 

(vi) growing to infinity of the Coulomb collision damping at a; — > a;^ + 0; 

(vii) non-harmonic composed waveform and signal dispersion and its smearing in cu; 

(viii) experimental requirements to electron distribution function to be Maxwellian and 
electrons to be coUisionless are intrinsically contradictory; 

(ix) geometrical (diffraction) effects (declination from the proposed one-dimensional prob- 
lem with plane waves since the length of the discharge tube is much more than its 
diameter) . 

The decrement of the amplitude damping at dispersion smearing of a wave can be 
calculated according to expression 



where /? = c or /3 = \Jv1 iov transversal or longitudinal waves, correspondingly; Acu is 
spectral width of the boundary exciting field £"(0,^). This smearing might be the main 
cause of the wave amplitude damping in coUisionless plasma both for longitudinal and 
transversal waves that imitates Landau damping. 

The return to Vlasov prescription of calculating relevant integrals in the principal value 
sense and the proper determination and using of mutually dependent boundary (or initial) 
conditions in the self-consistent manner allow to solve and remove all paradoxes of "Landau 
damping" . 

Thus, the right natural and simple but non-traditional ideology of waves in plasmas 
reduces to combine the method of Laplace transformation with the proper account for 
boundary/initial conditions, and Vlasov prescription of the principal value sense of relevant 
integrals. 

Mathematical correctness of the represented solution follows from the way of its con- 
struction. However one can analogously construct also other mathematically irreproachable 
solutions, for instance using "Landau functions" , that is with integrals in the principal value 
sense plus results of passing around the running with pi and p2 poles = —Pi/p2 in the 
lower or upper half-plane (in either case of fixed p° with two roots for the further calculated 
poles p^). 

It appears that the only criterion to pick out a physical solution from the set of math- 
ematically correct solutions to be only physical reasons, namely the definition of logarith- 
mically divergent integrals and related together boundary conditions in a such manner to 
avoid the paradoxes of divergent at x — > oo wave solutions, appearance of kinematical and 
backward waves, and waves not related to the boundary electrical field. The stabihty of 
the physical solution might be also considered as an additional requirement. 

Namely, these and only these criterions justify constructing the solution in terms of 
functions with integrals in the principal value sense as an optimal and apparently the 
unique variant. 

Violation of the uniqueness and existence theorem is related with singularity of the 
Laplace image of f^^^ due to logarithmic divergence of the integral in dvx at points Vx — 
— Pi/p2, where pi and p2 are the running Laplace transform parameters. 
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Attempts to establish functional constrains on the maximal value of the divergent inte- 
gral naturally prove to be artificial and^ can be justified only by special physical consider- 
ations. 

As a conclusion, our principally new approach in the plasma wave problem is the pre- 
scription of self-consistent calculation of the lacking boundary and initial conditions that 
gives a possibility to calculate the correct relative amplitudes of all relevant oscillatory 
modes and to eliminate the known paradoxes. Here acts a simple and robust principle: in 
the self-consistent problem there must be no perturbations which are not induced by (and 
as a consequence nonlinearly related with) the perturbing boundary electrical field. 

Transition to the real-valued boundary condition E{0,t) = Eocos^ut) in the asymptot- 
ical two- wave solution a ■ exp(ia;t — ikx) + a* ■ exp{—iujt + ikx) is equivalent to representing 
the solution as (cf. [2]) 

E{x,t) = EQCOs[Lut-kx + (p{k)] ; (3) 
fi{v,Vx,x,t) = F{v,Vx)cos[ujt- kx + tplk)] , (4) 

where E'q, (p{k), F{v,Vx), ipik) are real- valued amplitudes and phases which have to be 
determined. The substitution of E{x,t) and fi{v,Vx,x,t) into Vlasov equations at {p{k) — 
ip{k) = ±n/2 leads to the traditional equations 

^("."J = ±^f^l^. (5) 

(6) 
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with real roots ±k of Eq.® at u > ojl. If Eqs.© and (jH)) are satisfied, then E{x,t) (jSj) 
and fi{v,Vx,x,t) (jU are asymptotical solutions of Vlasov equations. 

The case of the totally imaginary k (Re {k) = 0, u < ul) in collisionless plasmas 
corresponds to the total wave reflection without energy dissipation. According to the afore- 
mentioned theory the solution reduces to a sum of terms, which equally exponentially damp 
at X — s> oo. As can be verified by direct substitution into Vlasov equations, elementary 
general asymptotical solution has the following form 

E{x, t) = EqC'^^^"" [cos{ujt + ^) + cos{ujt + ip)] , (7) 
fi{v,Vx,x,t) = e-^^^''[F{v,Vx)cos{ujt + if) + F{v,-Vx)cos{ut + i;)] , (8) 

m ovx oj^ -I- \kYv% 

with if — if) = ±7r/2 and dispersion equation being defined by Eq.ijH}. 

In dependence on the sense of improper integral in Eq. (jH)) one obtains different solutions 
k{(jj), with functions Q and (jH) being indeed the asymptotical solutions of Vlasov equations. 
From the real-valued character of the solution it follows directly the necessity to exclude in 
Eq.® additive imaginary constituents related to passing around the pole = uj/k, that 
is, to exclude the Landau prescription how to calculate this integral. From the physical 
symmetry of frequency and velocity of propagation of the forward and backward waves 
with respect to interchange k — /c it follows also that one has to calculate the integral 
in Eq.(jni) in the sense of principal value. This becomes more evident when one uses our 

^since real physical values are limited by different physical processes which are not accounted in the 
original equations in the given formulation of the problem. 
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approximate expression for the principal value of this integral. For calculation of the phase 
ip{k), amplitudes and initial and boundary values fi{v,Vx, x,0), fi{v,Vx,0,t) in order to 
take into account several modes ki^u), k2{uj), . . . one ought to use the general procedure of 
Laplace transform method with exclusion of kinematical and backward waves. 

Note, that solution at ^ uj/k tends to ±00, violating the initial perturbative 
condition |/i| ^ /o- However, one can assume from the physical point of view that there 
occurs some saturation of |/i| growth at periodically sign-alternating /i with limitation 
\fi{vx = oj/k ±£:)| < foivx = uj/k) at e — > 0, and presumably = uj/k) = with 

linear dependence /i on Vx in this region of unapplicability of the kinetic equation and 
possible discontinuity of dfi/dvx near /i ~ ±/o, so, at least, there are no physical reasons 
for appearance of any imaginary part in the integral (jH]) and in the function fi^ All 
these appear to be additional arguments in favor of calculation of the divergent integral 
in © in the principal value sense. The case of complex k in forward waves (Re A; > 
0, collision damping) was considered in our previous papers using the general Laplace 
transform method. In this case for real-valued boundary (and initial) conditions and u the 
solution is defined with a pair combination of exponents of the type a ■ exp(ia;t — ikx) + a* ■ 
exp{—iut + ik*x). 

Owing to the symmetry A; — > — in isotropic medium, the wave number k must enter 
into dispersion equation only through /c^. In coUisionless plasmas asymmetry of |/i| relative 
to the neighbourhood of the point Vx = uj/k±e leads to the possibility of appearing complex 
k, which define both damping and growing waves. Therefore this asymmetry leads to the 
paradoxical availability of waves (both in forward and backward directions) with exponential 
damping as well as with exponential growing due to the fact that asymptotic solution in 
this case contains both exponents ±i{ujt + kx) and ±i{ujt — kx). 

The fact of appearance of the indefinitely divergent integrals is the direct consequence of 
incomplete information in original differential equations. There must be additional exclu- 
sively physical (non- mathematical!) considerations at calculation of these integrals, since 
all different possible ways lead to different equally right (by their Laplace transform con- 
struction) mathematical solutions with correspondingly different dispersion relations k{(jj) 
and different wave amplitudes. But these solutions can be as physical as well as unphysical 
ones and have to be selected. 

To resume, the functions Q and (jH), as well as ((Tj) and (jH)), are evidently the very 
simple illustration of the existence of asymptotical non-damping or non-dissipative real- 
valued solutions for plasma wave equations. 

In the resonance region, where the condition /o(f , Vx, x, t) + fi {v, Vx,x,t) > is violated, 
the original, linearized or precise, equations are inapplicable. The divergence of /i at 
Vx = uj/k can not be removed when one uses expansion of the precise kinetic equation 
with the quadratic term in multiple overtones F„ cos[n{ujt — kx) + ipn] and En cos[n{ujt — 
kx) +ipn]- Thus, the original equations, linearized or not, are not sufficient to obtain single- 
valued solution of the problem without introducing into the kinetic equation (which details 
Liouville theorem) some additional physical specifications. 

In the theory of "Landau damping" PP velocity Vx can take, according to Cauchy the- 
orem, some arbitrary complex values at arbitrary deformation of the integration contour 
bypassing real- valued pole Vx = ui/k (with u > ul, Im (k) = 0) in the complex plane, that 
leads simultaneously to the appearance of exponentially damping and growing wave solu- 
tions with non- zero ±Im (k). According to |I] the latter is related with the small non- zero 

^In this case J \ fi\dv = \Ane\ can be treated as the real- valued concentration of perturbed electrons 
(|Ane| < n^). 
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values Im (k) resulting from the next derived quadratic dispersion equation. It evidences 
the intrinsic discrepancy and non-self-consistency of the Landau damping theory, since 
complex k, according to P], is a result of solving the dispersion equation, contrary to the 
before assumed real- valued k with the only one pole u/k on the real axis v^, thus the pole 
bypassing in the complex plane ought e.g. to be made now not along the half-circle (as 
it is presented in jTj), but along the total circle 0, or must not be made at all, with not 
clear grounds after all, for taking any value of k{uj). 

Conclusion 

The uniqueness theorem for solutions of the plasma wave equations, both for longitudinal 
and transversal waves at definite boundary and initial conditions, is violated due to the 
presence of logarithmically divergent integral in dispersion equations. To avoid the appear- 
ing of exponentially divergent solution terms, e. g. for considered plane waves in isotropic 
Maxwellian plasmas in the plasma half-infinite slab, one ought to treat this divergent inte- 
gral in the principal value sense. The physical requirements for the solution select the only 
solution out of all mathematically possible ones. We have presented elementary general 
real-valued solutions of Vlasov equations for collisionless plasmas at real-valued boundary 
and initial conditions in the form of a trigonometric function sum both for non-damping 
plasma waves (Im (k) = 0, u > uj^) and for the total reflection case (Re (k) = 0, uj < ujl) 
with dissipationless damping (evanescence) proportional to exp{—\k{uj)\x) . 

Finiteness of the wave solution in the neighbourhood of the point = uj/k ± e results 
by no means from energy dissipation in the sense of contradictory Landau damping but 
just from the purely kinematical effect of the electron density limitations with providing 
proper corrections to the kinetic equation. 

Existence of the general real-valued non-damping or dissipationless finite solutions of 
Vlasov equations is the direct proof of incorrectness of Landau damping theory. 
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